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The relativistic 2-body problem, much like the non- 
relativistic one, is reduced to describing the motion of an 
effective particle in an external field. The concept of a rela- 
tivistic reduced mass and effective particle energy, introduced 
some 30 years ago to compute relativistic corrections to the 
Balmer formula in quantum electrodynamics, is shown to 
work equally well for classical electromagnetic and gravita- 
tional interaction. The results for the gravitational 2-body 
problem have more than academic interest since they apply 
to the study of binary pulsars that provide precision tests 
for general relativity. They are compared with recent results 
derived by other methods. 

PACS number(s): 03.30.+P, 04.25.-g, 95.30.Sf, 97.80.Fk 



I. INTRODUCTION 

The notion of an effective relativistic particle describ- 
ing the relative motion of a 2-body system was first in- 
troduced in the context of the quasipotential approach 
to the quantum field theoretical eikonal approximation 
JIJ and bound state problem in quantum electrodynam- 
ics 0, Q. For a survey of subsequent developments - 
see B. Only later a general classical mechanical formu- 
lation of the relativistic 2-body problem was given || 
within Dirac's constraint Hamiltonian approach m (for 
a review - see 0). 

The central concept of a relativistic reduced mass is 
derived in this early work by observing that the total 
mass M of a 2-particle system should be substituted by 
its total centre-of-mass (CM) energy: 



w 1 
- = M + —E 

c c z 



M for 



Ma 



< 1 



(1) 



M = mi + m,2- 

This suggests using an energy dependent expression 
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M) , (2) 



for the relativistic generalization of the reduced mass /i. 
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Furthermore, if we determine the off-shell momentum 
square b 2 (w 2 ) for a pair of free particles as the solution 
for p 2 of the equation 



w 



= 1/ m\c 2 + p 2 



--b 2 (w 2 ): 



m\c 2 + p 2 => 

i — 2{m\+m 2 )c A w 2 + {m\ — m\) 2 c s ' 
Aw 2 c 2 



we find for the effective particle CM energy 



E w = c\J m^c 2 + b 2 (w 2 ) = 



w — mxc 



2w 



(3) 



(4) 



The interest in this notion of an effective particle was 
revived recently in |g] where a modified version of it was 
applied to the general relativistic 2-body dynamics and 
the relevance (and relative simplicity) of the dimension- 
less counterpart of the effective particle energy 



En 



2m 1 m 2 c 4 



(5) 



(which only makes sense for positive mass particles) was 
pointed out. On the other hand, the authors of || pre- 
ferred to work with the non-relativistic reduced mass /i 
rather than with the energy dependent quantity m w . 

The present paper was motivated by our wish to 
demonstrate the advantage of the original notion of rela- 
tivistic effective particle, cited in the beginning, in both 
the (classical) electromagnetic and gravitational 2-body 
problem. 

We begin in Sec. 2, by recalling the constraint Hamil- 
tonian approach to a relativistic particle system. A new 
justification is provided on the way for formula (^) for 
the energy dependent reduced mass identified as the co- 
efficient to the relative velocity in the expression for the 
effective particle 3-momentum in the CM frame. 

Sec. 3 is devoted to the electromagnetic interaction of 
two oppositely charged particles. 

The general relativistic gravitational two-body prob- 
lem (which continues to attract attention - see, e.g. M, 
PI, 0, |1, |H) is addressed in Sec.4. We compute 
the perihelion shift as well as the parameters of the last 
stable orbit in the gravitational case and compare with 
earlier results. 

In both cases we neglect the retardation effect: it is 
known not to contribute to the first post-Newtonian ap- 
proximation (see JTif) for which our results agree with 
previous calculations. The effects of the relativistic kine- 
matics, on the other hand, are computed exactly. 
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The possibility to take into account the finite velocity 
of propagation of interactions starting with the second 
post-Newtonian approximation within the effective par- 
ticle approach of this paper is discussed in the concluding 
Sec.5. 



II. VELOCITY SPACE FORMULATION OF THE 
CONSTRAINT HAMILTONIAN APPROACH TO 
THE RELATIVISTIC 2-BODY PROBLEM 

The mass-shell constraint for a free relativistic particle 
can be interpreted as a "Lorentz invariant Hamiltonian": 



Denote by P and w the CM momentum and the total 
energy: 



P 



■ Pi +P2, 



-P 2 c 



2„2 



(12) 



We shall exploit the fact that the difference (pi — ip 2 of 
the constraints ( |l0| ) is independent of the interaction i> 
to define the relative momentum 

P = M2Pi~MiP2, Mi+M2 = l (13) 

determining /Uj — /x 2 from the strong equation 

2Pp = 2<pi - 2(^2 = rn\c 2 + pf — in\c 2 - p\ 



1 

2A 



H = — {m 2 c 2 +p 2 ) »0, p 2 :=p 2 -p 2 Q . 



(6) 



Indeed, the equations of motion are obtained by taking 
Poisson brackets with H: 

i" = {x», H} = -p", p^ = {p", H} = Q (7) 
A 

for {x",pu} = 8S. 



Mi - M2 = 



2 4 



(14) 



(Thus, for unequal masses, the /i a depend on the Poincare 
invariant total energy square, so that the relation (13) is 
actually nonlinear.) The constraint ipi — tp% ~ together 
with (jlj) implies the orthogonality of p and P as an uni- 
versal kinematical constraint, readely solved in the CM 
frame: 



Here A is a Lagrange multiplier (assumed independent of 
x); it is linked to the choice of a time scale. The Hamil- 
tonian constraint gives rise to a singular Lagrangian 
through the Legendre transform, C = px — H for p deter- 
mined (as a function of x) from x — ^M- : 



C{x,x;\) (= px 



(8) 



Remark 1 It is important to remember - especially 
for the subsequent extension to a pseudo-Riemannian 
space-time - that the component p$ of the covariant 4~ 
momentum p^, canonically conjugate to the space-time 
coordinate x^ , is minus the energy, p = —E p /c (E p > 0; 
for a free particle E p /c = y/ m 2 c 2 + p 2 ). The equality 
p a = Ep/c is not generally covariant, it is accidentally 
valid for Cartesian coordinates in flat Minkowski space. 

For m 2 > A can be excluded from the condition 



x 2 -\ — ^-J=0 C(x,x) = -mc\[^x 2 . (9) 



d\ ~2 



Remark 2 For m = Eq. (jijj implies the constraint 
x and only the original expression (|; for the La- 
grangian remains meaningful. 

For a two particle system we introduce a pair of gen- 
eralized mass-shell constraints 



<Pa - g (P 2 a + m2 a + <p ) ~ ' a =1,2 
satisfying the (strong) compatibility condition 



{<Pl,¥>2} 



1>2 



d 
8x 2 



Pi 



d 

dxi 



$ = 0. 



(10) 



(11) 



(15) 



pP^O, (P=(w/c,0) => p=(0,p)). 



In order to solve the compatibility equation (|ll| ) we 
introduce the projection xj_ of the relative coordinate 
x i2 = xi — x 2 on the 3-space orthogonal to P: 



Xj_ = X12 



c 2 ^P. 



Its square, 



R 2 



12 



+ c 



{Px 12 f 



(16) 



(17) 



provides an invariant measure of the distance between 
the two particles in the CM frame. 

The general Poincare invariant solution of (|ll|) will be 
written in the form 



$ = $ (R,p 2 ,p R ; E w ) , where Rp R = px±. 



(18) 



The effective particle energy E w is singled out since, 
as we shall see, $ is a quadratic polynomial in E w in the 
cases of interest. One should also assume that $ — > for 
R — ► oo, thus making the separation of the mass terms 
in ( |l0| ) meaningful. 

The Hamiltonian constraint which replaces the equal- 
ity p 2 w b 2 (w 2 ) (jj^) in the presence of interaction is 
given by 



H := (fJ,2<Pi + M1V2) 



1 / 2 
2A> 



0. 



(19) 



Remark 3 The [i a defined in and have a 

simple expression in terms of the CM energies E a , 
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Ex = 



Eo = -- 



L 2 4 



2w 



2 4 



2w 



(20) 



a Galilean transformation that sends the CM velocity 
V = z& v, + ^ v, to zero gives vP M = vi - V = 



IT v > v 2 - 
CM frame we find 



v 2 — V = — tt£v.) In particular, in the 



and approach their non-relativistic values for 
M: 



E a 1 m( - m 

Ma = = o ± o 2 " C 



We no£e that there are precisely three ways to factorize 



b 2 into (f 
E w ft: ' 

2 ( W 2 ) 



El 



(^r + toc) /or i? egwa/ <o £a, i? 2 a^c? 



m 2 a c 2 = ^f-m 2 w c 2 , (a = 1,2). (22) 



They correspond to the i ( 2 ) subdivisions of the four 
zeros of the numerator of (0) into two pairs. 

This provides a fresh justification for the expressions 
(||) and (||) for the relativistic reduced mass and effective 
particle energy. We shall present yet another argument 
in favour of these expressions starting with a Lorentz co- 
variant concept of a relative velocity (cf. |D| Sec. 16). 

To simplify writing we choose units for which c = 1. 
The 4-velocities u a of the two particles and the CM ve- 
locity U are proportional to the corresponding momenta: 



Pa 



m a u a 



m a Uu a 



1,2; P = wU, 

~E a , U 2 = -l. 



(23) 



The CM energies E\ and E 2 (g0|) are thus related to the 
inner products of the velocities. 

We note that the constraints (|23| ) are equivalent to the 
kinematical constraint (|l5|). 

The 4-momentum of an effective particle with CM en- 
ergy e, u — (uq, u), Uq M = — e is introduced by defining 
first the relative 3-velocity u as follows. 

Let A = A(U) be the pure Lorentz transformation that 
carries the CM 4- velocity U into its rest frame. The con- 
ditions K^ U U V = <5p and positive definiteness determine 
the (symmetric) Lorentzian matrix A uniquely: 



A = 



U° 

-V 1 



U'Uj 

' T+TF' 



(24) 



We then find that the space parts of Aui and Ait 2 are 
proportional to the same 3-vector u: 



Aui 



m 2 . mi 
— u, Au 2 = u, 



w w 
which we shall identify with the relative 3-velocity: 



(25) 



miiii + miii2 



w/c 2 — *M 777 

' 1 ► 2* a=l,2. (21) 



E 2 

u = Ui 

m 2 



Eh_ 

mi 



-u 2 



m 2 



-Ui 



= 
w 
m i 



-u 2 . 



(27) 



The time component uq of the effective particle 4- 
momentum is determined by the condition: 



Uu 



.CM 



(28) 



Here e is given by (^|); for free particles it can be in- 
terpreted as the common value of the energy component 
—uio in the rest frame of u 2 and of —u 2 o in the rest 
frame of u\ which is a Lorentz invariant; in general, the 
following strong equation takes place 



UiU 2 ~ 



2m 2 



1) 



m 2 
2mi 



(29) 



-w 



2mim 2 



We identify, as usual, the space part of the effective 
particle 3-momentum p in the CM frame with the com- 
mon value of pf M and — p 2 M and define the proportion- 
ality coefficient between p and the relative 3-velocity u as 
the relativistic reduced mass m w : 



(Peff 



CM 
P 2 



m lU f M 



(30) 



where pf M can also be expressed in terms of the compo- 
nents of pi and pi in an arbitrary frame: 



P? M = A Pl = 



(E 2 +p 2 )p 1 -(E 1 +p° 1 )p 2 



w - 



with the non- 



(The relation (30) is consistent 
relativistic limit in which p CM = miVi = — m 2 v 2 = /itv.) 

The relativistic reduced mass m w defined by ( |30| ) is 
given, as anticipated, by (^|). The effective particle 4~ 
momentum Pej f is expressed in terms of the relative mo- 
mentum p (]l3|) and the CM 4-velocity U by 



Peff = E W U +p, 



n CM 
Peff 



(-E w , p), p 



CM 



(0, P). 



For positive mass particles it is convenient to write 
the Hamiltonian constraint (|l^) in terms of a dimen- 
sionless 4-momentum, corresponding to setting A = 
Am^c 4 (= \m 2 w ): 



u 



o 



u 



Uu 2 ) ui 



lit 







Uui) u 2 



1 + U° 



(26) 



H 



2\ 



(u 2 + $ 



l- e 2 ) 



0, 



(31) 



(Note that a similar procedure reproduces the non- 
relativistic relative velocity v = vi — v 2 . Indeed, 



<I> 



= $ (r, u 2 ,u r ;e) 



(32) 
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where the radial normalized momentum u r is dimension- 
less while r has dimension of an action: r = Rm w c, 
PR = m w cu r . 

The interaction function $ in ( fjll) is chosen to repro- 
duce the known interaction of a test particle in an exter- 
nal field. For the combined electromagnetic and gravita- 
tional interaction of two particles of charges e\ and e 2 we 
shall write: 

g = i-f 1+ ( 1 - 2 ik)„; + 4-(^n.o 

2A I V r J r r 2 1-2^ ) 

for r > 2a G (33) 

where J 2 is the square of the total angular momentum 
while a G (denoted in Q by a) is the gravitational cou- 
pling measured in units of action: 

CM 



J = r x u, r = m w cx± , ru r = ru, 



Here we have used angular momentum conservation 
which implies that tha effective particle moves in a plane 
orthogonal to J: 



r = r(cos v?, sin ip, 0) for J = (0, 0, J). 



(38) 



(34) 



We shall study the case of electromagnetic and grav- 
itational interaction (corresponding to a G — and to 
ha = 0, respectively) in Sec. 3 and 4 below solving in 
each case the resulting equations of motion. The result 
for the first relativistic ("post-Newtonian") approxima- 
tion agrees with (more complicated) traditional calcula- 
tions. Higher order corrections require taking into ac- 
count retardation effects which can be also done within 
the present 1-body approach, as discussed in Sec. 5. 
The Wheeler-Feynman non-local action for electrically 
charged particles Jl6| (see also 0) seems to provide a 
systematic treatment of retardation effects (to all orders) 
but its inclusion in the present framework is not obvious. 



This Hamiltonian constraint corresponds to the dimen- 
sionless interaction function 



$ = -2- 



r— 2a, 



-e 2 -2 



ha 



h 2 a 2 



G „,2 



2^u 



r — 2a G r(r — 2a G ) r 



eie 2 
ch 



(35) 



(a is positive for oppositely charged particles, studied 
in Sec. 3 below; it coincides with the fine structure con- 
stant , a -1 = 137.036, for e\ = — e 2 equal to the elec- 
tron charge). It is remarkable that the Poincare invariant 
constraints in flat phase space admit an interpretation in 
terms of an effective particle moving along geodesies in a 
Schwarzschild space-time. 

An advantage of our choice of A and of the variables 
r, u and e (instead of R = xJ M , p CM , and E w ) yield- 
ing the dimensionless Hamiltonian constraint ( J33| ) is the 
quadratic dependence of H in the (single !) energy pa- 
rameter e (instead of the two ^-dependent quantities E w 
and m w in (19). This allows to write down a Lagrangian 
for the (interacting) two-particle system using the stan- 
dard Legendre transform 

£(r, i, r, (p,;X) = ur — ei — H, ur = u r r + Jip, 



a r . \ ■ ha 



e = X 1-2^ i 



Ar 



, u,, , 

r / r 1—2- 



-, J=\r 2 <p (36) 



yielding 



A / f 2 



2 V 1 — 2^- V r J r 2A' 



■ ha 1 



or, varying in A and excluding it from the resulting con- 
straint, 

■ ha 



C = t A" 

r 



■ ha I ( „ „a„\ r 

t 1 - 2^ t 2 — 

r \\ r ) 1 - 2^ 



2,^2 



r ip 



(37) 



III. THE BOUND-STATE PROBLEM FOR TWO 
OPPOSITELY CHARGED RELATIVISTIC 
PARTICLES 

We start with the Hamiltonian constraint ([53]) in the 
absence of gravitational forces (i.e. with G = = a G ): 



I 2 

2XH = u 2 + — + l 



e+— I «0, 

r 



s 2 = ha 



eie 2 



(39) 



The canonical Poisson bracket relations {x£,pbu} — 
^abSfiiy imply the following non-zero brackets for the ra- 
dial and angular variables relevant for the planar motion: 



{r, u r } = 1 = {</>, J}. 



(40) 



The equations of motion derived from (39), (J4C|) read 
dH u r 



du r 



A 



-dJ=M^ J ^M-°- (41) 



The equation for the effective particle trajectory obtained 
by dividing r by <f> is independent of A: 



J 



J 



— — = u r = \\2ea, (l 



J 



(3 (42) 



where all variables - starting with — - are dimensionless, 
a j plays the role of a "classical fine structure constant" : 



J 



jeiegj 
cJ 



J 



(43) 



The bounded motion corresponds to the case when the 
expression under the square root has two positive zeros 
in — . This implies 



4 



< 13 = 1 -e 2 < a 2 , 
i.e., < -b 2 (w 2 ) < a^w 2 , 



(44) 



In fact, introducing the dimensionless inverse radius 
variable 



1 - ai J 



(45) 



we find 



^(^(-i^ + Wl- (46) 



The discriminant of the quadratic expression in y in 
the right hand side is positive whenever Eq. ( |4~4| ) takes 
place. (The positivity requirement for the two real roots 
in j implies that 1 — e 2 and 1 — a 2 are both positive.) 
We can then rewrite Eq. (j42|) in the form 



V(yi-y)(y-y2) v 



with j/i i2 = 1 ± e(0, J) 
where e(/3, J) plays the role of eccentricity: 



e 2 (/3,J) = l- 



1 



or. 







1-/3' 



0<e(e,J)<l, (48) 



the last inequality being valid in the domain (| 



J, i.e. 



for 



l<a 2 +e 2 , < ctj< 1, < e < 1 in the (a,, e) plane. 

Integrating Eq. ( j47| ) with initial condition y(<j) = 0) = 
yi (i.e., the orbit passes through the perihelion for cp = 0) 
we obtain 



y = 1 + e(/3, J)c OS (y'l-a 2 ^) . 



(49) 



To compare this result with the familiar non- 
relativistic elliptic orbit we first observe that for w given 
by (0) 6 @ * s ex P resse d in terms of the dimensionless 
measure e of the binding energy and the ratio v between 
the reduced and the total mass, 



E 
lie 2 

as follows 



M 



6=1 



£ (50) 



2 £ 



so that 

f3 = 1 - e 2 = -e (2 + (1/ + l)e) + 0(e 3 ). 



(51) 



(52) 



In the non-relativistic limit Eqs. (45) and ( |49| ) yield an 
elliptic trajectory 




COS( 



(53) 



The relativistic orbit fl49|), on the other hand, is not 
closed (except for e(/3, J) = 0). The perihelion shift 6<fr 
is given by 

<ty=27r((l -a 2 y 1/2 -l) = ira 2 j +0{a 4 j ) for a 2 < 1. (54) 
For a circular orbit we have 

e 2 (/3, J) =0 j3 = a], i.e., -e = f{a % n u) 



/(«» = - l- A /l-2i/(l- A /l-«», 



In general, for a bounded motion, we have the inequal- 



(47) ity 



0<-e</(« 2 ,^). 



(55) 



These results - and their derivation should be com- 
pared with the conventional approach that starts with 
the approximate Hamiltonian (see jL4|, Sec. 65, Problem 
2): 



H = 



2ju 



1 - 



l-Zi> p 2 
4 /i 2 c 2 



for p 2 = V \ 



I 2 

R 2 <* 



l + , p2+ ^ 
R V 2^ 2 c 2 



(56) 



& M \) 



A computation which involves a redefinition of i? ac- 
cording to the substitution R 2 => R' 2 = R^R — jjfj^ , 
and is certainly less transparent then the above, yields a 
trajectory whose parameters agree with ( [f9| ) up to (in- 
cluding) order a 2 (Velin G. Ivanov, Diploma Work, Sofia, 
2000). 



IV. GRAVITATIONAL 2-BODY PROBLEM 

The Hamiltonian constraint for the gravitational in- 
teraction of two (point) particles of arbitrary masses 
rni,m2, obtained from p3] ) for e\e2 = 0, can be in- 
terpreted as the condition that the effective particle 4- 
velocity u = (e, u) has unit mass in a Schwarzschild met- 
ric whose "radius" 2a G = m w cR w (of dimension of an 
action) is determined by the two masses: 



2A v 



0. 



(57) 



(A constraint of this type has been first written in |lq ] 
where, however, a more complicated metric was intro- 
duced, computed in a quantum field theoretic framework. 
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As observed in M the classical Schwarzschild metric gives 
a better approximation - in accord with our general pre- 
scription of Sec. 2.) Here the metric is expressed in terms 
of the radial variable r by 



5oo 



1 



1 



g' J UiUj = 1 - 2-± )K + — 



(<0), 



J 



The finite (bound state) motion belongs to the range 
2/2 < y < 2/i for which Pz{y) is non-negative. (The infi- 
nite interval y > 7/ , in which P?,(y) > as well, corre- 
sponds to falling on a centre.) The necessary and suffi- 
cient conditions for p and (3 for which all zeros of P3 are 
positive are the positivity of (3, 



(58) < p (= 1 - e 2 ) < 1, i.e., (m 2 + m\ <) < M 2 (65) 



where r, u r , J and a G are given by (|3J 

Proceeding to the Hamiltonian equations of motion we 
introduce the (J dependent) dimensionless coupling pa- 
rameter 



and the non-negativity of the discriminant: 



(66) 



R w m w c 2m\m2G 



J 



cJ 



(59) 



(i^ = 2^ being the energy dependent "Schwarzschild 
radius"). As we shall see shortly, p 2 < | for the bounded 
motion; by contrast, the counterpart a G /h of the elec- 
tromagnetic fine structure constant is rather big: for 
mi = 772 2 = Mq (the solar mass) it is of the order of 
10 76 .(The parameter p coincides with 2/j of |8).) 

The Hamiltonian constraint can be written in terms of 
p and - as 



H-± 1 
2A I 



J 



l-p- 



0. (60) 



The Poisson brackets (40) remain unchanged and we 
deduce as before the equations of motion 



. dH / J\ u r 



a/ = a^' (61) 



Introducing again a dimensionless variable propor- 
tional to the inverse radius (cf. (0)), 



J 

y = - 

r 



(62) 



we obtain the following (A-independent) differential equa- 
tion for the effective particle trajectory: 



We begin our discussion with the case of vanishing dis- 
criminant that gives rise to a circular orbit with 



2/1 



2/2 =: 2/c = 

3 o 
- + -p 3 
2 8 P 



1- y/r^y _ 

+ ^P 5 + <V), 



2/0 



- -2y c 
P 



obtained by solving the quadratic equation 

2 . -3/2 . 9 



(67) 
for 



P={py c -yc+p)yc 



27p 2 



(i-3 P y 



1 + 2^' 



A measure of the frequency on a circular orbit is 



dt 



= 4>/i = 



j 



i-p- 



J 



= ^{i- P y c ) (69) 



where we have used (J60|) , (|6l|) and (|62^ 

The Zasi (innermost) stable circular orbit, LSO (whose 
significance stems from the fact that gravitational radia- 
tion damping tends to circularize binary orbits j^] - see 
also H for a discussion and references) corresponds to 
the values 



= (! - py) u r = (py 3 -y 2 + py- 0) 

f3=l-e 2 . 



1/2 



(63) 



The energy independent coefficient p = p T will play 
the role of dimensionless expansion parameter (replacing 
the commonly used -). 

Eq. ( |63r) c an be solved in terms of Jacobi elliptic func- 
tions (cf. |19j] Sec. VII. 8). To begin with, we assume that 
all three zeros 2/0)2/i>2/2 of the cubic polynomial under 
the square root are positive reals 

^3(2/) ■■=py 3 -y 2 +py-f3=p(y-yo)(y-yi){y-y2), (64) 

o < 7/2 < 2/1 < yo- 



l 

7s 



—K = y c =yo> ft 



9' 



V8 
3 ' 



w 

— =M 



\ 



l-2v\l 




i *^2/3 1 



(70) 



for which both sides of the last inequality ( p6| ) vanish and 
all three zeros of P3 coincide. These values correspond 
to a limit point of local minima of what is called the 
"effective potential" V(r,J) that enters the expression 
for e 2 obtain from (60): 



6 



(71) 



We have 

± V ( r J)= J dP ^ 
dr ' r 2 dy 



(3py 2 -2y+p)=0. (72) 



It is the smaller zero of the second factor, y c (|67|), that 
corresponds to a minimum of V. The distance r c = — 
decreases when p — p , increases and attains its minimal 
value for the maximal possible value p* ( |70|) of p. 

The dimensionless binding energy e = — 
— {m w — p)jvm w of the LSO, evaluated from ( |5l| ) and 



M = 



= -I (l - v/l-2Kl-£*)) = 



■(1-6*) 1 + 1/ 



1-e* 



,(l-e* 



+ 0(^( e *) 3 ), (73) 



3 35 



The increase of |e 
Schwarzschild value |e 

, t-e* i 



— £*) (compared to its 
= 1 - e* ) by the factor 

1 + v^- + . . .J is coupled to a similar increase of the 

relativistic reduced mass 



1 - v\e* 
for ^- = M(l+i/e*). 



M(l + ^e*) 

if" 
c 



(74) 



The expansion parameter x = (a G o;) 2 ^ 3 of |l3| ] takes 
for the LSO its Schwarzschild value x* — | - sec (70). 



Returning to the "true radius" R = — — (measured in 
units of length) of LSO we find, in accord with that 
it is smaller than its Schwarzschild value R* 



R*=K^=K(^l-Rs(i--[i-^))- (75) 



Accordingly the angular frequency uj* for the relativis- 
tic two-body LSO is bigger than its Schwarzschild value 



for < p 2 < -. 

3 

(The coefficient to p 8 is chosen in a such a way that for 
p 2 = \ inequalities ^% become equalities.) We shall 
compute the zeros of P3 as expansions in p 2 taking into 
account the fact that /3 is, according to (|77|), (at most) 
of order p 2 . 

For the largest root one finds the following expansion 

Vo = - - P (1 + P 2 - P + 2p 4 - i P 2 p) + 0(p 7 ). (78) 
P 

The two smaller roots are then computed from the re- 
lations 



2/1 + 3/2 = Vo 

P 

p(l + p 2 -/3 + 2p 4 -3p 2 /3)+0(p 7 ), 
pyoym^P- (79) 

The result is 
t/1,2 = f (1 + P 2 - (3 + 2p 4 - 3p 2 /3 ± A) + 0(p 7 ), 



A > 0, (80) 



2 x 2 

M = 

p 2 -4/3+2 ( o 4 -6 j o 2 /3+5/7 6 -16 j o 4 /3+5 / 9 2 /3 2 + 0( i o 8 ). (81) 

The solution of Eq. ( |63|) satisfying y(0) = Wi is expressed 
in terms of the elliptic sine function (see EQ]): 



2/W - 2/2 2 /V AT \ ^ ; 

= sn \K - V P\yo ~ 2/2) 77, K 

2/1-2/2 V 2 



(82) 



The module square, k 2 , of the elliptic functions is ex- 
pressed as a ratio of differences of roots of P 3 : 



k 2 = 2/l"2/2 



Ap 2 



2/0-2/2 l-f/^l+^-^ + l^A 



+ 0(p 8 )= (83) 



V ( 1 • %r • 4 P 4 - |p 2 /3 



( P 2 -4/3 + 5p 4 -18p 2 /3)+0(p 8 ); 



Mc 2 



1 + ve* 



(76) 



We now proceed to the general case in which the rela- 
tion (|6§) between f3 and p 2 becomes an inequality 



< < 



Tip 2 



[l-3p 2 f 2 -l 
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< 



1 2 

4 P 



45 o 1 , 



(77) 



4K(k 2 ) is the rea/ period of sn(x,K) and cn(x,k): 
sn(K,k) = l,cn(K,k) = 0, 



7o v /(l-^ 2 )(l-^ 2 ) 2 V 2 2 



(84) 



We now proceed to computing the perihelion shift 5(f) 
of the effective particle. The change A<j> of <j> for a full 
turn and the shift 6cj> are given by 
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4iT 



2/2 \/P3{y) vw - 2/2) 



= 2tt + <5</>. (85) 



This is an exact formula. For small p the elliptic mod- 
ule k is also small (according to j83|)) and we can approx- 
imate using (p4), by 



2tt [ 1 + -p 



4K = 2;t I 1 
A , / : *» 



4 
,2 



64 



1 



0(k 6 ) 

Combining this with the expansion 

(p(Vo - y2)y 1/2 = 



(86) 



l-^ 2 (lV-/?) + ^A 



\p 2 (l + |(3p 2 -2/3) 



-1/2 



v 

4 



^o(p 6 ; 



we end up with an expression for 
powers of A cancel out: 



in which the odd 



AK 



Vp(vo - 2/2) 



- 2tt = 



3tt 



P 2 ( 1 + ^( 7 ^ 2 - 4 ^) + ^ 6 ) 



To compare with earlier calculations fllO| ], | pT[ , p^ ] 
one again uses the expansion ( |52j) of /3 in terms of the 
dimensionless binding energy (pCj). The result clearly 
agrees with the first post-Newtonian approximation. The 
missing fourth order (in p) term ^j-vp 2 (/3 — |p 2 ) can be 
shown to correspond to retardation effects. 



V. CONCLUDING REMARKS 

Formulae for particle trajectories in a relativistic 2- 
body system (including recoil effects) have been derived 
with the same ease as for a test particle problem in a 
Coulomb or Schwarzschild potential. The expression e = 
1 



e + |e (|3lJ) for the CM energy per unit mass e 



in terms of the dimensionless measure e 



lie- 



of 



the binding energy (e < for finite motion) has been 
deduced as a straightforward consequence of the relation 
(H) between the CM energy, the total mass M = rrii+rri2 
and E: 



w 
Mc 2 



1 



JL 

M 



TO1TO2 



(mi + m-iY 



(89) 



(This simple and natural derivation should be compared 
with the rather involved argument of Sec. 4 of ||] yielding 
the same result. It thus provides a posteriori justification 
of our definition of a relativistic effective particle - on top 
of the arguments presented in Sec. 2.) 



A systematic way to compute higher order corrections 
has been worked out in the quantum case [§|J|], and can, 
in principle, be applied to the classical limit as well. It 
is, however, desirable to have a consistent classical algo- 
rithm for calculating retardation effects. 

It is known (and follows by comparing results of the 
previous sections with earlier calculations) that these first 
contribute to order e 2 (i.e. to a* in the electromagnetic 
case or to p A (~ 4r) in the gravitational case). It turns 
out that it is quite feasible to include such corrections 
within the effective 1-body approach developed here. We 
shall indicate how to do this by modifying the effective 
particle Lagrangian. 

It follows from Eq. (26.23) of || (or from Eq. (65.7) 
of [Q) that the retardation effect in electrodynamics (to 
order Jy) is accounted for by multiplying the interaction 

term (i.e., — ) by the velocity dependent expression 



1 - ~ ( viv 2 + i(vir)(v 2 r) 



(90) 



(where we are using our dimensionless velocities v a , cor- 
responding to -v a , a = 1, 2, in the above references). 

Similarly, for the gravitational case (following Sec. 103 
of |l5|]), one has to multiply the interaction term (i.e., 
— £-) in the Newtonian Lagrangian with the velocity de- 
pendent expression 



k — k 

G EM 



\ ( 3 (vi -v 2 ) s 



(91) 



We recall that a simple minded application of the above 
procedure would give spurious first order effects that 
should be eliminated. 

Noting that the CM velocities of the two particles are 
expressed in terms of the relative velocity v by 



E w 

Vi = — V, V 2 : 

Ei 
we find 



E u 
~Ei 



-v for v„ 



1 E 2 



E a 



v=^u (92) 



2EiE 2 

. wE w 



(93) 



E1E2 



V 



In the above approximation one should replace the en- 
ergy dependent factors with their non-relativistic limits 



El 



E1E2 



wE w 
EiE 2 



1 for e —i- 0. 



(94) 



We end up with the following simple expressions depend- 
ing on the effective particle 3-velocity v and position r 
only: 



k EM = !+ 2 U ( V2 + V r 



(95) 



8 



It is clear that these new terms will account for the dif- 
ferences between our result ( |B8l) for 5cj> (obtained neglect- 
ing retardation effects, i.e., under assumption k G = 1) 
and the result in [H (in order p 4 , p 2 /3). It is, how- 
ever, challenging to develop the corresponding calcula- 
tional scheme in order to be able to compare the third 
post-Newtonian approximation tackled in Ea]. 
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